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This paper reviews the work of Wu and Lin on the Union Jack lattice Ising model. This model 
is of interest as it one of the few to display re-entrant phase transitions. Specifically we re-examine 
their result for the general anisotropic sublattice magnetisations, comparing these with the works of 
Vaks, Larkin and Ovchinnikov, and our own numerical simulations. We discuss the disagreements 
found in both sublattice predictions including non-zero antiferromagnetic results and a rotational 
variance. We will then suggest additional conditions and modified formulae that will allow valid 
results to be produced. 

PACS numbers: 64.60 De, 64.60 Bd, 64.70 qd 



I. INTRODUCTION 

The Ising model is the prototypical model of phase 
transitions and as such is greatly studied [l|. Is is made 
up of sites interconnected along a lattice of "bonds". 
Each of these sites can have spin of value ±1. The Union 
Jack lattice is obtained by adding alternate diagonals to 
a square lattice as shown in Figure [TJ This lattice is made 
up of two sublattices, the cr-lattice with sites with eight 
intersite "bonds" (filled circles) and the r-sublattice with 
sites with four of these "bonds" (hollow circles). 




FIG. 1: The Union Jack Lattice. (From [2|) 

The Union Jack lattice Ising model is of particular 
interest as it is one of the few exactly solvable models 
which exhibits a re-entrant phase transition |3(. A so- 
lution was presented for the general anisotropic model 
on this lattice for the c-sublattice in Q and then later 
for the complete lattice in The solution showed that 



the Union Jack lattice Ising model is equivalent to the 
free fermion model. The spontaneous magnetisation for 
the free fermion model had been computed in [5|. The 
Union Jack lattice Ising model has been greatly studied 
since these solutions, with the work being extended in 
d, 0] to study a mixed spin lattice. 

This paper focuses on the numerical simulations of the 
work in the thesis on the general anisotropic Union Jack 
Ising model in [8]. To start, in Section [TT] we will state 
the required background for the paper. In this section we 
will also present the sublattice prediction functions from 
the work of @, 0] ■ In Section IIIII we will briefly discuss 
the method used in our simulation program, along with 
the qualitative modelling of the Union Jack Ising model. 
We will then go onto present a theoretical analysis of 
the work presented in 0, 0] and @ and compare these 
results with those from our numerical simulations. We 
will show that the prediction functions of @, 0] do not 
accurately model the results of many of the anisotropic 
systems with invalid results being produced. Specifically 
we will identify the cases where the r-sublattice predic- 
tions are physically implausible and those which are not 
consistent under rotation of the lattice. In the conclu- 
sion in Section ITVl we will state the additional conditions 
required to allow the predictions to produce valid results 
and the causes of the rotational variance. 



II. THE UNION JACK LATTICE 

In this section we will briefly state the equations for 
spontaneous magnetisation for the general anisotropic 
Ising model on the Union Jack lattice as presented by 
Wu and Lin 0, H|. For convenience we adopt the no- 
tation of Wu and Lin Q and label the nearest interac- 
tion strengths as — J r , which are defined to be one of six 
values: — J±, — J2, —J3, —J a, —J, — J'- The resultant 
Boltzmann factors are given by 
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where (3 — 1/ksT, kg is the Boltzmann constant and 
a,b,c and d are the four sites surrounding a r site. This 
equation produces sixteen possible factors, which can be 
reduced by symmetry to eight distinct expressions 0: 



In it is shown that the Union Jack lattice is equivalent 
to an eight-lattice models with weights given by (fTJ) . This 
eight lattice model satisfies the free fermion condition |9j . 
The spontaneous magnetisation of a free fermion model 
was given by Baxter in Q. As such, the spontaneous 
magnetisation for the cr-sublattice is 
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The calculation for F + and F_ is a little more involved. 
We start by calculating 
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We can relate F + and F_ with the following formula from 
[Hi, allowing us to get values for each variable, 
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We can compute the overall nearest neighbour magneti- 
sation by taking the mean of the two sublattice magneti- 
sations 

M = h{a) + {r)). 



71 = —Ul + U>2 + W 3 + LU4 

7 2 = uji — U2 + UJ3 + UJ4 

7 3 = uii + lo 2 — L03 + uja 

-f4 = Ull + L0 2 + LO3 - UJ4. 

The critical point (s) of this system are 

n 2 = 1 

or equivalently, 

uji + oj 2 + uj 3 + LU4 = 2 max{wi,w 2 , w 3 , w 4 } . (4) 

The r-sublattice magnetisation is given by [4( as 

(r) = {a) [A 1234 (K){F+ + F_) + 4 2341 (iT) (F + - *"-)] • 

(5) 

We can see that this equation is a multiple of the cr- 
sublattice value. In §5§ from j4j], 
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and 



G-{(3J) = cosh2(/3J! + (3J 3 ) + cosh2(/3J 2 - /3J 4 ). 



A. Classification of phases 

In 0, a classification of the phase of the cr-sublattice 
is presented. It is based on the following energy values: 

-E x = J + J' + I Ji + J2 + J3 + Ja\ 
-En = -J-J'+ |Ji - J 2 + J 3 - J 4 | 
-JSs = - J + J' + \Ji -J2-J3 + Ja\ 
-E A = J - J' + I Ji + J 2 - J 3 - J 4 | . 

The sublattice is in: 

i) a ferromagnetic phase when Ei < E2, E 3 , E4; 

ii) a antiferromagnetic phase when E2 < Ei, E 3 , E4; 

iii) a metamagnetic phase when E 3 < E%, E2, E4 or 
E4 < Ei, E2, E 3 . 

As the temperature rises, depending on the relative 
strengths of the interactions J r , the occurrence of a phase 
change is signified by one or more of the equations in (U) 
being realised. A re-entrant phase transition occurs if 
any one equation admits two solutions. 



III. RESULTS 

In our work in Q we presented a theoretical and nu- 
merical analysis of the results derived in Q • In our theo- 
retical analysis we compared these results against those of 
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3] to identify systems to further investigate numerically. 
The results of both the theoretical and numerical analysis 
will be presented below. In [|| we also presented analysis 
of a mean-field approximation of the Union Jack lattice 
using two approaches. The first used a set partially un- 
coupled predictor equations, both functions of (a). The 
second approximation used coupled predictor equations 
that are functions of (a) and (r). Qualitatively we saw 
that our Mean Field models showed good correlation with 
the isotropic ferromagnetic systems. In anisotropic an- 
tiferromagnetic systems as well as anisotropic ferromag- 
netic systems with a re-entrant phase transition, the cor- 
relation was poor between the Mean Field results and the 
theoretical predictions. 

As the theoretical results are for an infinite lattice, we 
used a Monte Carlo Markov Chain method with periodic 
boundary conditions to numerically simulate the systems. 
Our Monte Carlo algorithm was the Metropolis-Hastings 
algorithm [l|, [Io| . For our simulations we will apply this 
algorithm to a lattice of 100 sites by 100 sites. This lat- 
tice size has been chosen as it is small enough to have a 
reasonable run time while being large enough to suppress 
the finite size effects. The computer program was cali- 
brated against exact results for the general anisotropic 
triangular lattice for both the average magnetisation [ll| 
and the three-site correlator [T2|. With this calibration 
we have confidence in the accuracy of the simulation re- 
sults. The code for this simulation program can be found 
in @. 

To classify the phase of the er-sublattice we used a al- 
ternative approach to that given in Section III Al By us- 
ing the individual 7 terms from ([3]) we note that the 
sublattice is in a ferromagnetic phase when 71 < 0, a 
antiferromagnetic phase when 72 < 0, a metamagnetic 
phase when 73 < or 74 < and a disordered phase 
when 71727374 > 0. The critical temperature(s) can also 
be determined when any of the 7 functions move from 
being negative to being positive, or vice versa. 

A. Isotropic ferromagnetic 

To start with we look at the isotropic ferromagnetic 
system. Here the interactions for our system will be J = 
J' = J n = 100/c b- In our initial plot of this system the r- 
sublattice prediction was a factor of two higher than 
expected. Upon analysis of the equation we found that 
due to the symmetric interactions of the system A1234 — 
^2341 = 1, F + F- = and (r) = 2F + (a). As such we 
adapted ([5]) to the following form 

, . , . [A 1234 (K)(F + + F_) + A 2341 (K)(F + - F_)} 

(v = w 2 ■ 

(7) 

As this is a minor adjustment we will continue to refer 
to the result as that of Wu and Lin [|| . The plot of our 
numerical simulation results against the adapted predic- 
tion results of Wu and Lin [2, 13] are shown in Figure [5] 
Intuitively this is the type of graph we would expect for 
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FIG. 2: Plot of simulation results for an isotropic 
ferromagnetic system on the Union Jack lattice. Here 
the theoretical predictions are shown with the lines and 
the numerical results are shown with the points. 

this type of system. As expected there is strong agree- 
ment between the two theories with the phase transition 
and critical temperature being the same. There is a high 
correlation between our simulation data and both sublat- 
tice prediction functions. There is some noise around the 
critical temperature, though it is of a small magnitude 
when compared to the other results. After the noise part 
of the data, around 400-600 Kelvin the simulation results 
again follow the prediction with a high correlation. 

B. Anisotropic metamagnetic 

Next we move on to look at an anisotropic meta- 
magnetic system where J n — lOks, J — lOOks and 
J' = — lOOfcs. The graph we obtain when we plot our 
simulation results against the predictions of Wu and Lin 
0,3 is shown in Figure[3]below. As can be seen from the 
graph, the prediction functions of 0, H| show a non-zero 
magnetisation on the cr-sublattice In this system the 74 
is negative at low temperatures Q. The equation for 74 
is 

74 = UJl + uj 2 + W3 - W4, 

and we note from examination of this equation that 

LJi = lu (+ H ) is the dominant term. From the spin 

configuration, this represents we would expect an antifer- 
romagnetic phase with the average magnetisation for the 
CT-sublattice being zero. When we look at our simulation 
results for this sublattice we see that they have average 
spin zero which follows more closely the results we would 
expect. 
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FIG. 3: Plot of simulation results for an anisotropic 
metamagnetic system on the Union Jack lattice. Note 
that the simulation results do not follow the curves of 
the prediction functions. 



FIG. 4: Plot of simulation results for an anisotropic 
antiferromagnetic system on the Union Jack lattice. 
Note that the simulation results do not follow the 
curves of the theoretical predictions which are 
physically impossible. 



C. Anisotropic antiferromagnetic 

Next we study an anisotropic antiferromagnetic sys- 
tem. The chosen system of this type will have horizontal 
and vertical interactions of J n = lOOfcs and diagonal in- 
teractions of J = J' = — 100k b- The simulation results 
are plotted against Wu and Lin's 0, 13] predictions in Fig- 
ure [4] below. Using similar analysis as in Section MI Bl 

we find that U2 = U)(-\ 1 — ) is the dominant term for 

the magnetisation of the c-sublattice. This configuration 
would agree with the result of [3j in that it predicts an 
ordered antiferromagnetic phase and the critical temper- 
atures agree. However as we can see in Figure [H the 
prediction function ([3]) produces a non-zero magnetisa- 
tion for this phase where we would expect zero magneti- 
sation. 

When we examine the r-sublattice prediction func- 
tion ||7J) we notice that it produces physically implau- 
sible results. That is, the predicted magnetisation of 
the r-sublattice is greater than the possible spin val- 
ues, ±1, for this Ising model. As in the isotropic fer- 
romagnetic case we have symmetric interactions leading 
to A1234 = ^2341 = 1 and F+ F- = 0. As in this case 
F + > 2, we can see that this is the cause of the implau- 
sible results. 

The correlation between the simulation data and the 
predictions of @, E| is again poor. We see that our simu- 
lation results are centred around the zero magnetisation 
level, as we would intuitively expect. However it can be 
seen that the simulation results for the c-sublattice do 
not show a difference between an ordered phase and a 
disordered phase. 



D. Anisotropic ferromagnetic 

We now move to examining anisotropic systems with 
the cr-sublattice classified as ferromagnetic at low tem- 
peratures. This type of system is particularly interesting 
as it contains systems with staggered interactions and 
those with re-entrant phase transitions. We will also dis- 
cuss rotational variance. First let us present a system 
containing staggered interactions, with the interactions 
of the system shown in Figure [5] defined as J n — — lOOfce 
and J = J' = 92k b- As we can see from the prediction 
functions of @, 0] the average magnetisation have oppo- 
site sign. This agrees with the theoretical predictions for 
the critical temperature and state from [3j. It can also 
be noted that the overall magnetisation for the complete 
lattice has magnitude zero until around 200 Kelvin. The 
numerical simulation results have a high correlation with 
the theoretical predictions and follow the curves closely. 
There is noise after the critical temperature on both sub- 
lattices but this quickly reduces down to a small level for 
higher temperatures. 

Having covered staggered interactions we now go onto 
looking at systems that contain a re-entrant phase tran- 
sition. The system we studied was one with horizontal 
and vertical interactions of J n — lOOfcs and diagonal in- 
teractions of J = J' = —92ks- The graph of this system 
is shown in Figure H>] This system is predicted by [H 
to start in a ferromagnetic phase, move to a disordered 
phase, then into an ordered antiferromagnetic phase and 
then finally into a disordered phase again. As we see 
from the prediction functions of [2j, |4j , initially the theo- 
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FIG. 5: Plot of simulation results for an anisotropic 
ferromagnetic system on the Union Jack lattice with 
equal horizontal and vertical interactions. 
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FIG. 6: Plot of simulation results for an anisotropic 
ferromagnetic system on the Union Jack lattice with 
equal horizontal and vertical interactions. 



ries agree with a ferromagnetic phase being shown. The 
critical temperatures predicted by both theories agree in 
all three values. With further analysis we see that ini- 
tially 71 is the dominant term, leading to a ferromagnetic 
phase. From the second critical temperature to the third 
we see that 72 is now the dominant term. However as we 
saw in Section 1111 CI while we would expect a zero mag- 
netisation, there is instead a non-zero value for both sub- 
lattices, with the r-sublattice having value greater than 



one. 

In our numerical simulations we see that there is good 
correlation up to the second critical temperature but the 
simulation does not show the re-entrant phase transition. 
As we have discussed in Section IIII CI both the ordered 
antiferromagnetic phase and disordered phase would have 
average magnetisation zero, and so are indistinguishable 
from this point of view. 

Looking at more a general anisotropic system, we 
see that there exist systems which have rotational vari- 
ance. As an example of these systems we first look at 
the system with horizontal interactions of J\ — J3 = 
lOO/cs/0.9 2 , vertical interactions of J 2 = J4 = 100£;b/0.9 
and diagonal interactions of J = J' = lOOfcs, and 
compare with the system with horizontal interactions of 
Ji = J3 = lOOfcs/0.9, vertical interactions of J2 = J a = 
100/cb/0.9 2 (i.e. the lattice rotated through 90 degrees). 
The graphs of our simulation results are shown in Figure 
[7J below. As we see when we compare the simulations 
results, the data forms similar curves and has a simi- 
lar phase transition at equal critical temperatures. In 
comparison to the predictions of Wu and Lin 0, 3 > we 
see that at higher temperatures the simulation results 
follow all three curves with good correlation. At lower 
temperatures, below about 200 Kelvin, we see that the 
a prediction still has good correlation for both systems. 
This suggests that the r-sublattice prediction is incor- 
rect. Our simulation results show that rotation of the 
lattice should not have an effect on the results of the 
system. When we analyse (J7J) we see that 

p p _ COsh 2 2/? ./Horizontal ~ CQsh 2 2/3J Vort ical 
COSh 2 2/3 ./Horizontal + COSh 2 2/3 ./vertical 
^4-1234 = 1 

COSh 2/3 ^Vertical ,„s 

A 234i — — r do T ■ \°) 

COSh i/JJHorizontal 

Both the rotational variance and disagreement between 
the r-sublattice and the simulation results are a result of 
©. 

A further result can be seen if we now take a system 
similar to the previous example but with negative diag- 
onal interactions. For an example of this type of sys- 
tem, we will look at the system with horizontal inter- 
actions of Ji = J3 = lOOfcs/0.9 2 , vertical interactions 
of J2 = J4 = lOO/cs/0.9 and diagonal interactions of 
J = J' = — 100k b, and the same system rotated through 
90 degrees. The graphs of our simulation results are 
shown in Figure [5] Again we see as in Figure [JJ that the 
simulation results are very similar to each other. How- 
ever, in this case we see that below 10 Kelvin all three 
simulation results are lower than the predictions. Af- 
ter 10 Kelvin the simulation results again move up to 
the prediction curves, following the a prediction in both 
cases. We see that in both cases the simulation results 
and predicted results show the same critical tempera- 
ture and phase transition. This oddity maybe due to the 
simulation being performed on a finite system, while the 
theoretical predictions are for an infinite lattice. 
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FIG. 7: Plot of simulation results for an anisotropic ferromagnetic system on the Union Jack lattice with horizontal 
and vertical interactions that are not equal and positive diagonals, (a) shows the system with interactions 
Ji = Ji = lOOfcs/0.9 2 , J2 = J4 = lOOfcs/0.9 J = J' = 100/cb. (b) shows the rotated system with interactions 
Ji = J 3 = lOOfcs/0.9, J 2 = J 4 = 100/cb/0.9 2 J = J' = 100k B - 




FIG. 8: Plot of simulation results for an anisotropic ferromagnetic system on the Union Jack lattice with horizontal 
and vertical interactions that are not equal and negative values for the diagonal interactions, (a) shows the system 
with interactions J\ = J3 = IOO/cb/0.9 2 , J% = J4 = lOOfcs/0.9 J = J' = — 100/cb- (b) shows the rotated system with 
interactions Ji = J 3 = 100fc B /0.9, J 2 = J 4 = 100/c B /0.9 2 J = J' = -lOOfcu. 
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IV. CONCLUSION 



We have seen that the re-entrant phase transitions of 
the Union Jack Ising model can not be seen when only 
considering the average magnetisation. This is due to the 
transition being from unordered to an ordered antiferro- 
magnetic phase which both have an average magnetisa- 
tion that is identically zero. In addition, the prediction 
of the cr-sublattice given in Q requires additional condi- 
tions to have agreement with our numerical simulations. 
It is possible to classify the phases of the system by exam- 
ining the 7i terms of equation ([3]). The predictions with 
the current conditions produce non-zero magnetisations 
for non-ferromagnetic systems. However, if we impose 
the conditions that we only use the prediction formula 
^) when 71 < or 71 72 73 74 > and it is zero outside 
those conditions, the results now work for all systems. 



The prediction for the r-sublattice, given in has 
additional issues in the general anisotropic lattice. Ini- 
tially the prediction given in [J] is too large by a factor 
of two, but can be easily corrected as shown in (|7J). The 
cr-sublattice prediction is a factor of the r-sublattice pre- 
diction and so a zero magnetisation would be observed 
for the antiferromagnetic and metamagnetic phases if the 
above conditions are applied. The rotational variance 
seen in Figure[7]and Figure[S]can be eliminated by rewrit- 



ing ([6]) as follows 
A 123 4(K) = 



sinh2(/3J 1 +/3J 3 ) 



A 23 4i(K) = 



y/2G\-{(3J) sinh 2/3 Ji sinh 2/3 J 3 

sinh2(/3J 2 +/3J 4 ) 
v/2G 2 - 03 J) sinh 2/3 J 2 sinh 2/3 J 4 ' 



where 



Gi-(/3J) = cosh2(/3Ji + pj 3 ) +cosh2(^J 2 -/3J 4 ), 
G 2 _(/3J) = cosh2(/3J 2 + /3 Ji) + cosh2(^Ji - /3J 3 ) 

This change also stops the disagreement at low tempera- 
tures in these systems. Applying both these changes and 
the conditions for the cr-sublattice gives agreement with 
the simulation results. 

Further work is required to understand the implica- 
tions of the conditions suggested above on the work of 
Wu and Lin In this paper they use a similar ap- 
proach to the one presented here for the checker-board 
lattice, so an investigation into this lattice would be use- 
ful to see if similar results are obtained. Also as many 
papers, such as @, 0], extend the results presented here. 
These should also be studied for similar inconsistencies. 
Finally, we observed a disagreement between the theo- 
retical predictions and simulation results at low temper- 
atures in our analysis of the system presented in Figure 
[8] This disagreement remains after the application of 
the above conditions. As such, further investigation into 
these systems is required. 
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